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Abstract 

Assume that there is a set of monic polynomials P„(z) satisfying the second-order difference equation 

A(s)P„(z(s + 1)) + B(s)P n (z(s)) + C{s)P n {z{s - 1)) = A„P„0(s)), n = 0, 1, 2, . . . , N 

where z(s),A(s),B(s),C(s) are some functions of the discrete argument s and N may be either finite or 
infinite. The irreducibility condition A(s — l)C(s) is assumed for all admissible values of s. In the finite 
case we assume that there are N + I distinct grid points z(s), s = 0, 1, ...,N such that z(i) ^ z(j), i ^ j. 
If N — co we assume that the grid z(s) has infinitely many different values for different values of s. In 
both finite and infinite cases we assume also that the problem is non-degenerate, i.e. A„ 7^ A m , n 7^ m. 
Then we show that necessarily: (i) the grid z(s) is at most quadratic or q-quadratic in s; (ii) corresponding 
polynomials P n (z) are at most the Askey- Wilson polynomials corresponding to the grid z(s). This result 
can be considered as generalizing of the Bochner theorem (characterizing the ordinary classical polynomials) 
to generic case of arbitrary difference operator on arbitrary grids. 
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1. Introduction 



General orthogonal polynomials (OP) P n (x) can be characterized by the 3- term recurrence relation [6] 

P„+\{x) + b n P n (x) + u n P n (x) = xP n (x) (1.1) 

with initial conditions P a = 1, Pi = x — b 

The polynomials P n (x) are monic polynomials, i.e. P n {x) = x n + 0(x n ^ 1 ). 

It is well known [1] that all polynomials solutions P n (x) of the second-order differential equation 

<j{x)PH{x) + t(x)P^(x) = X n P n (x) (1.2) 

are classical orthogonal polynomials (COP), i.e. Jacobi, Laguerre, Hermite and Bessel polynomials. In (1.2) it 
appears that a(x) and t(x) are polynomials such that deg(a(x) < 2, deg(r(x) = 1. This result is known as the 
Bochner theorem [5]. 

It is natural to consider generalization of the Bochner theorem replacing the second-order differential operator 
with the second-order difference operator. Namely we are seeking polynomial solutions P n {z) of the problem 

A(s)P n (z(s + 1)) + B(s)P n (z(s)) + C(s)P n (z(s - 1)) = X n P n (z(s)), n = 0, 1, 2, . . . , N (1.3) 

where z(s),A(s),B(s),C(s) are some functions of the discrete argument s and N may be either finite or infinite. 
The irreducibility condition A(s — l)C(s) ^ is assumed for all admissible values of s. In the finite case we 
assume that there are N + 1 distinct grid point z(s), s = 0, 1, . . . , TV such that z(i) ^ z(j), i ^ j. If N = oo we 
assume that the grid z(s) has infinitely many different values for different values of s. In both finite and infinite 
cases we assume also that the problem is non-degenerate, i.e. A„ ^ A m , n ^ m. We assume also that there are 
polynomial solutions of all degrees n = 0, 1, . . . , N (i.e. we assume that the polynomial P n (x) always has exact 
degree n for all n — 0, 1, . . . , N.) 

Askey and Wilson [2] discovered orthogonal polynomials (the Askey- Wilson polynomials, or briefly, AWP) 
which satisfy equation (1.3) for quadratic z(s) = as 2 + bs + c or q-quadratic grid z(s) = aq s + bq~ s + c, where q is 
some parameter such that |g| ^ 1. Finite-dimensional case (i.e. when there exists only N mutually orthogonal 
polynomials n = 0, 1, . . . , N — 1) corresponds to the so-called q-Racah polynomials [9]. 

In [7] it was shown that the only OP satisfying (1.3) for AW-grids are the AWP. Leonard showed [10] that in 
the finite-dimensional case the only OP satisfying (1.3) are the q-Racah polynomials. For further development 
of the Leonard result and its new algebraic interpretation see, e.g. [17]. In [8] Ismail obtained more strong 
result : he showed that all polynomial (i.e. not necessarily orthogonal, ab initio) solutions of the equation (1.3) 
for the AW-grid are AWP. In the finite-dimensional case Terwilliger obtained the result that the AW-grid is the 
most general for polynomials satisfying (1.3). 

So far, the open problem was: in the infinite-dimensional case characterize all possible grids z(s) for which 
polynomial solutions of the equaiton (1.3) are obtained. In this paper we solve this problem and show that there 
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are no grids more general than AW-grids. Hence all polynomial solutions for (1.3) should be orthogonal Askey- 
Wilson polynomials. Although for the finite-dimensional case the problem was effectively solved by Terwilliger 
in [18], we present here the finite-dimensional version of the generalized Bochner theorem as well. The main 
reason is that our method of proof is essentially different and deals directly with difference equation (1.3) for 
polynomials whereas in the Terwilliger paper [18] another (a purely algebraic) approach is presented. 

2. Finite-dimensional case 

In this section we show that if N is finite then the problem is essentially equivalent to the Leonard theorem 
[10]. 

Indeed, consider (N + 1) x (N + 1) tri-diagonal matrix J which acts on a basis e k , k = 0, 1, . . . , N by 

Je k = C(k + l)e fe+ i + B(k)e k + A(k - l)e fe _i (2.1) 

It is assumed that C(N + 1) = A(—l) = which means merely that that the matrix J acts in linear space of 
dimension N + 1. We will assume the nondegeneracy condition: 

C(i)A(i - 1) ^ 0, i = l,2...JV (2.2) 

Find the eigenvectors v^ k \k = 0, 1, . . . , N of the matrix J, i.e. 

J W « = \ kV M 

with some eigenvalues Afe. We assume that all eigenvalues are distinct: Aj ^ Xj if i ^ j. Then all vectors 
i/ fe ' , k = 0, 1, . . . , N are independent and we have 

N 
s=0 

where v ks , s = 0, 1, . . . , N are components of the vector in the basis e s . For them we have relation 

A(s)v kiS+ i + B(s)v ks + C(s)w fe!;s _i = \ k v ks . (2.4) 

Now we can identify components v ks with P k (z s ), i.e. we merely put v ks — P k (z s ) for all values k, s — 0, 1, . . . , N. 
Then difference equation (1.3) coincides with (2.4). 
Consider transposed Jacobi matrix J* defined as 

J*e k = A{k)e k+1 + B(k)e k + C(fc)e fc _i (2.5) 

and corresponding eigenvalue vectors v*^: 

J*v*W = \ kV < k \ fc = 0,l,2,...,JV (2.6) 
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Vectors v*^ can be expanded in terms of the same basis e s 



N 



v*W=J2v* ks e s (2.7) 

s=0 

From elementary linear algebra it is known that in nondegencrated case (i.e. if ^ Xj for i ^= j) the vectors 
v k and v*W arc biorthogonal: 

N 

(v k ,v*^) = Y / Vk s v* s = if k^j (2.8) 



s=0 

Introduce now the diagonal matrix M wich acts on basis e s as 

Me s = ft e s , s = 0,1, 2,..., N (2.9) 

where 

"° = 1 (2 - 10) 

Note that all are well defined due to nondegeneracy condition (2.2). 
It is elementary verified that 

J* = M~ X JM, (2.11) 

and hence 

v *(k) =M"V fe) , fc = 0,l,...,iV (2.12) 

(inverse matrix M _1 exists due to nondegeneracy condition (2.2)). Relation (2.12) allows one to rewrite biorthog- 
onality condition (2.8) in the form 

N 

w s v ks v js = 0, if k^j, (2.13) 



s=0 

where 

^v^ni (2.i4) 

In terms of polynomials P n (x) this relation becomes 

N 

Y J w s Pk{z{s))P j {z{s)) = Q, if k^j (2.15) 

But relation (2.15) means that P n {x) are polynomials which are orthogonal on a finite distinct set of points 
z(s),s = 0,1,..., N with discrete weights w s ^ 0. By general elementary theorems concerning orthogonal 
polynomials [6] this means that polynomials P n (x) should satisfy a three- term recurrence relation 

P n+1 (x)+b n P n (x)+u n P n - 1 (x) =xP n (x), n = 0,l,...,N (2.16) 

The roots x s , s = 0, 1, . . . , N of the polynomial Pn + i(x) coincide with spectral points: 

z(s) = x s , s = 0,l,...,N 
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Thus we proved that (under some nondegeneration conditions) polynomials P n (x) satisfying relation (1.3) on 
a grid z(s) for finite N arc orthogonal with respect to discrete weight function (2.15) and satisfy three-term 
recurrence relation (2.16). 

Now we are ready to relate our results with Leonard's approach to dual orthogonal polynomials [10]. 

Recall relation between nondcgenerated Jacobi matrices and orthogonal polynomials (see, e.g. [6]). Let 
K be an arbitrary Jacobi matrix of dimension N + 1 x N + 1. In some finite-dimensional basis d n it can be 
presented as 

Kd n = a n d n+1 + (3 n d n + j n d n -i (2.17) 
with some (complex) coefficients with nondegeneracy property 

N 

I] 7*^-1^0 (2.18) 

i=l 

Construct eigenvectors ir^ of the matrix K: 

Kir^=z k ir (k \ k = 0,l,...,N (2.19) 

We assume that all spectral points Zk, k = 0, 1, . . . , N are distinct: z k ^ Zj for k ^ j. Expand eigenvectors tt^ 
in terms of basis d n : 

N 

7T (fc) =^2w ks d s 
s=0 

with some coefficients 7Tfc s . For these coefficients we have from (2.19) the recurrence relation 

7 s+ l7Tfe iS+ l +/3 S -K ks + a s _l7Tfe iS _l = Z k TTks, k, s = 0,1, . . . , N (2.20) 

It is assumed that a_i = -Jn+i = 0. Then, for each value k, starting from nko we can find recursively all further 
coefficients nki,7Tk,2, ■ ■ ■ , ^k,N- We can always normalize 7Tfeo = 1, k = 0,1, . . . ,N. Then it is clear from (2.20) 
that TTks is a polynomial of degree s in argument Zk- 

Introduce polynomials T n (x) satisfying three- term recurrence relation 

ln+iT n+ i{x) + f3 n T n (x) + a n _iT„_i(a;) = xT n (x) (2.21) 

with initial conditions a_iP_i = 0, Po(x) = 1. Then relation (2.21) defines n-degree polynomials T n (x) = 
n n x n + 0(x n ^ 1 ) with the leading coefficient 

1 

7i72 • ■ • 7« 

(this leading coefficient is well defined and nonzero do to nondegeneracy condition (2.18)). From general theory 
of orthogonal polynomials it follows that polynomials T n (x) are orthogonal on a finite set of points Xk [6] 

N 

^2^kT n (x k )T m (x k ) = 0, n^m (2.22) 

k=0 
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where x k are roots of the polynomial Tn+i{x). 

We can thus associate orthogonal polynomials T n (x) with expansion coefficients of eigenvectors of the arbi- 
trary nondegenerated Jacobi matrix K: 

T s (z k )=ir ks (2.23) 

Return to our polynomials P n (x) satisfying relation (1.3). We showed that these polynomials are orthogonal 
and correspond to the Jacobi matrix K whose matrix coefficients can be restored from recurrence relation 
(2.16): 7„ = 1, (3 n = b n , a n = u n+ \. On the other hand, we have the Jacobi matrix J defined by (2.1). By just 
described recipe, we can associate with this Jacobi matrix corresponding orthogonal polynomials Y n {x). These 
polynomials satisfy 3-term recurrence relation 

A(n)Y n+1 (x) + B(n)Y n (x) + C(n)Y n ^(x) = xY n (x) (2.24) 

Now it is seen that polynomials P n (x) and Y n (x) are related as 

P n {z{s)) = Y s (X n ) (2.25) 

We thus have a duality property coinciding with that introduced by Leonard [10]: there are two systems of finite 
orthogonal polynomials and two finite sequences z(s) and A n such relation (2.25) holds. Our nondegeneracy 
conditions: all z{s) and A„ are distinct and matrices J, T are nondegenerated coincide with similar conditions 
in the Leonard paper. Hence we can conclude 

Theorem 1 Under nondegeneracy conditions the finite- dimensional case of relation (1.3) generates at most 
finite Askey-Wilson orthogonal polynomials (Racah and q-Racah polynomials in other terms). 



3. Infinite-dimensional case. Reducing to a more simple problem 

In this section we start to analyze the infinite-dimensional case. We first derive some restrictions upon the 
coefficients A(s) 7 B(s),C(s). 

In what follows we will assume that polynomial solutions P n (z) of the equation (1.3) are monic, i.e. P n (z) = 
z n + 0(z n ^ 1 ). This is not restriction of our problem, because it is possible to divide all terms in equation (1.3) 
by a (nonzero) leading coefficient of the polynomial P n (z). 

First of all we observe that eigenvalues A„ can be shifted by an arbitrary constant A„ — > A„ + const. Such 
shift leads to adding a constant to the coefficient B(s). Using this observation we always can choose A„ in such 
a way that 

A =0 (3.1) 

In what follows we will assume that condition (3.1) is fulfilled. We will also assume that the eigenvalue problem 
(1.3) is non-degenerate, i.e. 

A„ ^ A m , n ^ m (3.2) 
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The grid z(s) is also assumed to be non-degenerate, i.e. 

z(s 1 ) ^ z(s 2 ), si ^ s 2 



(3.3) 



Parameter s takes infinite number of integer values: s — s , s + l, s + 2, . . . where so is either finite or s = — oo. 
In the first case we deal with semi- infinite grid z s , whereas in the second case we have the grid which is infinite 
in both directions. 

Taking the case n — in (1.3) we see that A(s) + B(s) + C(s) = 0. Hence we can rewrite equation (1.3) in 
the form 

A(s)AP n (z(s)) - C(s)VP n (z(s)) = X n P n (z(s)), (3.4) 
where we use the standard notation [13] 

AF(s) = F(s + 1) - F(s), VF(s) = F(s) - F(s - 1) 

for any function F(s) of the argument s. 

Assume that polynomials P n (z) have the expansion 

P n (z) = Z n + J2Uz i 
i=0 

with some coefficients £„j. Then for n = 1 we get from (3.4) 

A(s)Az(s) - C(s)Wz(a) = XiQi(z(s)), 
where Qi(z) = z + £io = Pi{z). By induction, it can be easily shown that 

A{s)Az n {s)-C{s)Wz n {s) = \ n Q n {z{s)), n = 0,l,2,..., (3.5) 

where Q n {z) is a monic polynomial of degree n. 

Vice versa, assume that property (3.5) holds for some z(s), A(s), C(s) with Q n (x) being a set of monic 
polynomials in x of degree n. Then there exists a set of monic polynomials P n {x) satisfying equation (3.4). 
This statement is almost obvious and follows from the observation that on the given grid z{s) and for any monic 
n-th degree polynomial T n (x) the expression A(s)AT n (z(s)) — C(s)\7T n (z(s)) is again a n-th degree polynomial 
in the argument z(s) with the leading coefficient A„. Hence, it is possible to choose a polynomial P n {x) with 
the property (3.4). 

Consider now condition (3.5) for n — > n + 1: 

A( s )Az n+1 ( s )-C( S )Vz n+1 ( s ) = X n+1 Q n+1 (z( s )), n = 0,l,... (3.6) 
Multiplying (3.5) by z(s) and subtracting (3.6) we get another set of conditions 

A 1 (s)z n {s + l) + C 1 (s)z n {s-l) = R n+1 (z{s)), n = 0,l,... (3.7) 
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where A\(s) — A(s)Az(s), C\(s) = C(s)Vz(s), The polynomials R n (z) are n-th degree polynomials R n (z) = 
uj n z n + 0(z n ^ 1 ), where iv n = A„ — \ n -i- Note that due to non-degeneracy condition (3.2) we have u„^0 and 
hence every polynomial R n {z) has exact degree n. 

Consider first two conditions (3.7) corresponding to n — and n = 1. These two conditions can be considered 
as equations for two unknowns Ai(s), C\(s). Solving these equations we have 

Ra(z(a)) - z(a - 1)11^(8) 



M*) 



z(s + 1) - z(s - 1) 



n f , R2(z(s))-z(s + l)R 1 (z(s) 

C1{S) = z(s + l)-z(s-l) (3 ' 8) 

Note that these expressions are well defined for all possible s because, by non-degeneracy condition, z(s + 1) ^ 
z(s-l). 

Hence conditions (3.7) can be rewritten as 

R2(z(s))Y n - R x (z{a))z(a - l)z(a + l)y„_i = R n+1 (z(s)), n = 2, 3, . . . , (3.9) 

where 

_ z"(, + l)-z"Q,-l) 
Yn ~ z{s + \)-z{s-\) ■ ( ' 

Introduce the variables 

u = z(s - l)z(s+ 1), v = z(s-l) + z(s + l) 

Clearly Y n is a symmetric polynomial with respect to z(s — 1), z(s + 1) and hence it can be expressed in terms 
of variables u, v only. Indeed, it is easily verified that Y n satisfy the recurrence relation 

Y n+1 = vY n - uY n _ u Y =0,Y 1 = 1. (3.11) 

This allows us to find an explicit expression for every Y n in terms of u, v. For example, Y2 = v, Y3 = v 2 — u 1 Y^ = 
v 3 — 2uv etc. 

Return to condition (3.5). We have explicit expressions for coefficients A(s), C(s): 

Riz(s - 1) - R 2 



A(s) 



(z(s + 1) - z(s - l))(z(s + 1) - z(s)) : 



r( , ) = Riz(s + 1)-R 2 

1 > ( Z (s + l)-z(s-l))(z(s)-z(s-l)) { ■ > 

Hence we have 

n—l n—1 

A(s)Az(s) - C(s)Vz(s) - Yl * n ~ k ~\z(s + - ^1^-1 - RiYk) = - ^ ^+i^(s)""' s " 1 

k=0 fe=0 

(in the last equality we have used (3.9)). It is seen that this expression is indeed a polynomial of degree n with 
non-zero leading coefficient A„. Thus conditions (3.5) and (3.7) are equivalent and we can use only more simple 
condition (3.7) for further analysis. 
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4. Functional equation for the grid z(s) 

From (3.9) and (3.11) we find the conditions 

Rn +2 (z(s)) =vR n+1 (z(s)) -uRnizis)), n = 2, 3, . . . (4.1) 

These conditions form a system of linear equations for two unknowns u, v. Consider the first two equations 
corresponding to n — 2 and n = 3. There are two possibilities: 

(i) these equations arc not independent. Then we should have Ri+\{x) = r(x)Ri(x), i = 1,2,3 where t{x) 
is a linear function. By induction, we then have R n {x) — i?i(a;)T™ _1 (x), n = 1,2,... for all n, where both t(x) 
and Ri(x) are linear functions in x. Now from (4.1) we have the condition 

t 2 (z(s)) -vt(z{s)) +u = (4.2) 

or, equivalently, 

(t(z(s)) - z(s + l))(r(z( S )) - z(s - 1)) - 0. (4.3) 

From (4.3) and (3.12) we see that in this case either A(s) = or C(s) = for every admissible s. But this 
contradicts our non-degeneracy assumption A(s — l)C(s) ^ 0. Thus the case (i) should be excluded from 
consideration. 

(ii) these equations are independent. Putting n = 2,3 in (4.1) we obtain a linear system of equations for 
unknowns u, v from which we find 

7T 6 (z(s)) ^6(Z(S)) 

where TTi(x) are polynomials of degrees < i: 

7Tq = i?3 — R2R4, 7T7 = R4R3 — R2R5, 7Tg = R3R5 — R\ 

Thus u, v are some rational functions in the variable z(s). In what follows we will sometimes replace the grid 
z(s) with independent variable x (this is possible because the grid z(s) takes infinitely many different values). 

We first prove an important statement concerning possible solutions of the system of non-linear difference 
equations of the form 

z{s-l) + z{s + l)=T 1 (z(s)), z{s-l)z(s + l) = T 2 (z{s)) (4.5) 
where Ti_ 2 {x) are some rational functions. 

Lemma 1 Assume that the system (4-5) has a solution z(s), s — s n ,s + 1, . . . with infinitely many non- 
coinciding values z(si) ^ z(s 2 ) if si 7^ s 2 . Then there are two possibilities: 
(i) either 

m , \ a 2 x 2 +a4,x + a 5 a 3 x 2 + a 5 x + a 6 

T i{x) = h—, ; . T 2 (x) = 5- ; 4.6 

aix z + a 2 x + a 3 a.\x/ + a 2 x + a 3 

with some constants a>i, i = 1, . . . , 6. 
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In this case variables z(s),z(s + 1) satisfy equation 

$(*(«), *(s + l)) = 0, (4.7) 
where $(x,y) is a non-reducible symmetric bi-quadratic polynomial: 

$(x, y) = aix 2 y 2 + a 2 xy(x + y) + a 3 (x 2 + y 2 ) + a A xy + a 5 (x + y) + a 6 (4.8) 

or 
(ii) 

a 2 x + ot A a 3 x + a 4 „ /x a 2 x + a A a 3 x + a 4 

±\{x) = , l 2 (x) = (4.9) 

aix + a 3 a.\x + a 2 ct\x + a 3 a\x + a 2 

with some constants on,i — 1, . . . , 4 such that a 2 ^ a 3 . In this case variables z(s), z(s + 1) satisfy equation 

a 1 z(s)z(s + 1) + a 2 z{s) + a 3 z(s + 1) + a A = (4.10) 

Remark. The case (ii) formally corresponds to a special case of (i) when polynomial &(x,y) can be decom- 
posed as a product of two polynomials of the first degree in both variables x, y. 

Proof. Obviously, system (4.5) is equivalent to the statement that both z(s + 1) and z(s — 1) are roots of 
the quadratic equation 

A 2 (z(s))z 2 ±1 + A 1 (z(s))z s±1 + A (z(s)) - 0, (4.11) 

where Ai(x) are non-zero polynomials having no common factors. 
Introduce two polynomials in two variables: 

Wx{x,y) - A 2 (x)y 2 + Ai(x)y + A (x), W 2 (x,y) = Wi(y,a;) - A 2 (y)x 2 + A^x + A (y) 

Equations W\{x,y) = and W 2 (x,y) = define two algebraic curves in complex variables x,y. From (4.11) it 

is clear that both curves contain infinitely many common distinct points (x n , y n ), n = 1, 2, By the Bezout 

theorem this is possible only if these curves either coincide or have a common component. 
The polynomial W\{x,y) has degree 2 in variable y and hence there are two possibilities: 

(i) Wi(x,y) is irreducible, i.e. it cannot be decomposed into irreducible polynomials of a lesser degree in y. 

(ii) W\{x,y) can be presented as a product of two polynomials, each of degree 1 in variable y: W\{x,y) — 
(ei(x)y + e 2 (x))(e 3 (x)y + e 4 (x)) with some polynomials Rei(x),i = 1, ... ,4. 

We consider these two possibilities separately. In the case (i) we have that the polynomials W\{x,y) and 
W 2 (x,y) are both irreducible. Hence, by the Bezout theorem, they should coincide: 

W 1 (x,y) = W 2 (x,y) = W 1 (y,x) (4.12) 

But condition (4.12) means that the polynomial W\{x, y) is symmetric in variables x, y. This is possible only if 
all polynomials Ai(x), i = 0, 1, 2 have degree < 2 in variable x. Hence, the most general expression for W\{x, y) 
in this case is symmetric bi-quadratic polynomial in x, y: 

Wi(x,y) = aix 2 y 2 + a 2 xy(x + y) + a 3 (x 2 + y 2 ) + a 4 xy + a 5 (x + y) + a 6 (4-13) 
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with some constants a^, i = 1, . . . , 6. We thus have 

A 2 (x) = a\x 2 + a 2 x + a 3 , A\(x) = a 2 x 2 + a^x + a 5 , A (x) = a 3 x 2 + a 5 x + a e 



and 



a 2 x 2 + a^x + «5 . . . . . . . a 3 x 2 + a$x + 



Ti(x) = -A^/A^x) = - ' - ° , T 2 (x) = A (x)M 2 (x) 



o^x 2 + + a 3 ' aix 2 + a 2 x + a 3 

giving expression (4.6). 

Consider now the case when polynomials Wi(x, y), W 2 (x 1 y) have a nontrivial common component which 
doesn't coincide with both these polynomials. Clearly, this is possible only if Wi(x,y) can be decomposed into 
two polynomials linear in y: 

Wi{x,y) = ( ai (x)y + b 1 (x))(a 2 (x)y + b 2 (x)) (4.14) 
with some polynomials ai t2 (x),bi t2 (x). By definition W 2 (x,y) — Wi(y,x) and hence we have also 

W 2 {x,y) = (ai(y)a: + &i(y))(o2(y)a: + 62(»)) (4.15) 

Without loss of generality we can assume that a\ (x)y + bi (x) is a common component of two curves W\ (x, y) = 
and W 2 (x, y) = 0. Comparing (4.14) and (4.15), we can conclude that there are two possibilities: 

(i) cither a 1 (x)y + b x (x) = a 1 (y)x + bi(y); 

(ii) or a 1 (x)y + h(x) = a 2 (y)x + b 2 (y) 

In case (i) we have that variables x, y satisfy symmetric polynomial relation 

aixy + a 2 (x + y) + a 4 = (4.16) 

with some constants ai,a 2 ,ot4. Substituting x = z(s),y = z(s + 1) into (4.16) we find from (4.16) that there 
are only two non-coinciding points z(so) and z(s + 1). For all further points we find that z(s + 2j) = z(sq) 
and z(s + 2j + 1) = z(s + 1) for all integer j. But this contradicts our assumption that there arc infinitely 
many distinct points belonging to the curves. Thus the case (i) is impossible. 

In the case (ii) the polynomials ai :2 (x), b\ i2 (x) should be linear in x and we have that z(s), z(s + 1) satisfy 
the relation 

a 1 z(s)z(s + 1) + a 2 z(s) + a 3 z(s + 1) + a 4 (4.17) 

where a 3 ^ a 2 in order to prevent impossible case (i). This case corresponds to (4.9) and (4.10). Thus the 
Lemma is proven. 

It is interesting to find explicit solutions in both cases (i) and (ii) of the Lemma. The case (i) corresponds 
to a parametrization of symmetric Euler-Baxter bi-quadratic curve $(z(s), z(s + 1)) = with &(x, y) given by 
(4.8). This problem was already solved by Baxter [4] in his famous solution of the 8- vertex model. Explicitly 

z(s) = K0(/3 lS + /3 o ) (4.18) 
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with some parameters k, /?i,/?o- Here <j>(z) is an even elliptic function of the second order (i.e. having exactly 
two poles in the fundamental parallelogram). Recall that (up to an arbitrary factor) any even elliptic function 
of the second order can be presented in the form [19] 

m = I" = Pr + i (4-19) 

Recently it was shown that the elliptic grid z(s) described by (4.18) appears naturally in theory of biorthogonal 
rational functions with the duality property [15], [16]. Note that in a special case a\ — a 2 = the rational 
functions Ti(x),T 2 (x) become linear and quadratic polynomials. In this case solution for z(s) is expressed in 
terms of elementary functions (see below) . 

For the case (ii) of the Lemma the solution can be easily found in terms of elementary functions of s (we 
will not describe these solutions in details because they can be obtained from the elliptic solutions by a limiting 
procedure). 

Now return to condition (4.1) and consider first the case (i) of the Lemma. We can rewrite (4.1) in the form 

1 = v(x)R n+1 (x)/R n+2 (x) - u(x)R n (x)/R n+2 (x), (4.20) 

where 

. a 2 x 2 + a A x + a 5 a 3 x 2 + a 5 x + a 6 

v(x) = , u(x) = 5 

aix z + a 2 x + «3 a\x z + a 2 x + 013 

Assume first that a\ 7^ 0. Then for x — > 00 it is seen that rhs of (4.20) tends to which contradicts to lhs of 

(4.20). Thus necessarily a\ — 0. Assume now that a.\ = and a 2 ^ 0. Then again for x — > 00 we obtain from 

(4.20) the condition (recall that R n {x) = w n s" + 0(x n ^ 1 ), where u n = A„ — A„_i) 

1 = -Wn+l/Wn+2 

whence A n = \ n +2 for all n = 2,3,.... But this contradicts our condition of non-degeneracy of the spectrum 
A„. We thus have necessarily a\ = a 2 = 0. But in this case v(x),u(x) become polynomials of the first and 
second degrees: 

v(x) = —£,x — 77, u(x) = x 2 + rjx + ( (4-21) 
where £ = 0:4 /as, r/ — a 5 /as, ( — a 6 /a 3 and equations for the grid become 

z(s - 1) + z(s + 1) = -&(s) - rj, z(s-l)z(s + l) = z 2 (s)+ V z(s) + ( (4.22) 

with arbitrary complex parameters £,r], (. Equivalently, variables z(s),z(s + 1) belong to a non-degenerating 
conic (i.e. ellipsis, hyperbola or parabola): 

z 2 {s + 1) + z 2 {s) + r](z{s + 1) + z(s)) + £z(s)z(s + 1) + C = (4.23) 

which is symmetric with respect to z(s), z(s + l) (this means that the plot of this conic in Cartesian co-ordinates 
x = z(s), y = z(s + 1) is symmetric with respect to the line y = x). Equations (4.22) and (4.23) were studied in 



12 



[13], [11], [12]. In these works it was shown that all nondegenerate solutions of these equations can be presented 
in the form 

z(a) = C iq s + C 2 q- S +C (4.24) 

or 

z(s) = C 2 s 2 + C lS + C Q (4.25) 

or 

z(s) = (-l) s (C 2 .s 2 + da + C Q ) (4.26) 

with some constants Co, C\, d- The first case (4.24) occurs if £ = q + where g ^ ±1 (i.e. £ ^ ±2). The 
second case (4.25) occurs if £ = — 2 and the third case (4.26) occurs if £ = 2. All these cases exhaust possible 
types of the Askey- Wilson grids [12]. 

Note that when Cid = in (4.24) we obtain so-called exponential grids, say z(s) = dq s + Co- Similarly, 
when d — in (4.25) or (4.26) we obtain the linear grid: z(s) = C\s + Cq or z(s) = (— l) s (Cis + Co). However, 
in these case the conic (4.23) becomes degenerated - it divided into two lines. This corresponds to the case (ii) 
of the Lemma (see below). 

Now substituting v(x),u(x) into (4.1) we obtain that for arbitrary given polynomials Ri(x), R 2 (x) one can 
construct uniquely the polynomial R n {x) — uj n x n + 0(x n ^ 1 ) with leading coefficient satisfying the recurrence 
relation (which easily follows from (4.1) for given v(x)u(x)): 

uJ n +2 + + ui n = 0, n = 2, 3, . . . (4.27) 

General solution for u„ can be easily found from (4.27): if £ = q + q 1 ^ ±2 we have 

uj n = dq n + G 2 q- n (4.28) 

with arbitrary d, G 2 - If £ = — 2 then 

w„ = Gm + G (4.29) 

and if £ = 2 then 

w„ = (-l)"(Gm + G ) (4.30) 

From w n = \ n — \ n -i we can easily reconstruct the spectrum A„ which has the same functional dependence on 
n as the Askey- Wilson grid z(s) has on s. 

Finally, we should consider the case (ii) of the Lemma. In this case similar considerations lead to conclusion 
that «i = 0. Then condition (4.10) becomes 

ct 2 z(s) + a 3 z(s + 1) + a 4 = (4.31) 

which describes exponential or linear grids z(s). Thus the case (ii) can be considered as a degeneration of 
the case (i). It should be noted that linear and exponential grids are described by non-symmetric recurrence 
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relations (4.31) (with respect to z(s),z(s + 1)). The reason is that in this case the conic (4.23) is degenerated 
to a two straight lines each of which is non-symmetric . 

We see that in both cases (i) and (ii) of the Lemma solutions R n {x) of the recurrence relation (4.1) are 
indeed polynomials of exact degree n. Hence, by previous considerations, we obtain a unique set of polynomials 
P n (x) which are solutions of equation (3.4). From explicit form (3.12) of the coefficients A(s),C(s) we can 
conclude that they coincide with those defining the Askey- Wilson polynomials [13], [11], [12]. 

We thus proved that under some non-degeneracy conditions, the only admissible grid is the Askey- Wilson 
grid and corresponding polynomials P n {x) coincide with the Askey- Wilson polynomials. 



5. Concluding remarks 

The authors of [13] exploited relations (3.4) as a starting point in their approach to construction of the Askey- 
Wilson polynomials. In a slightly different manner, Magnus in [11], [12] derived relations (4.23) from the 
following requirement: find all the grids z(s) and y(s), such that for any polynomial P n (x) of degree n we have 
the property 

P n (z(s + 1))-P n (z(s)) 

Z(S+1)-Z(S) -^-liV^h (5.1) 

where T„_i(x) is a polynomial of degree n — 1 and relation (5.1) should be valid for all n = 1,2, .. . and for 
infinitely many distinct values s of the grids z(s) and y(s). Relation (5.1) can be also presented in the form 

V s P n (x)=T n _ 1 (y(s)) (5.2) 

where T> s stands for "discrete derivation" operator which acts on the space of function f{x) as 

' K ' z(s + l) - z(s) 

For the AW-grid it was known that the operator D s satisfies property (5.2). Magnus proved that these grids 
are the only preserving property (5.2). 

On the other hand, it was noted in [13] that if polynomials P n {x) satisfy the AW-cquation (3.4) then the new 
polynomials T n (x) obtained from P n (x) by (5.1) also satisfy AW-equation (3.4) but with different coefficients 
A(s), C(s). This property can be considered as a covariance of the Askey- Wilson equation (3.4) with respect to 
the discrete Darboux transformation (see, e.g. [14]). 

However the property (5.1) cannot be directly derived from equation (3.4) if the grid z(s) is not concretized. 
This is why derivation of the necessity of the AW-grid for equation (3.4) is not quite elementary and needs 
rather involved technique which was demonstrated in the present paper. 
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